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PACS 36 . 20 . Ey - Conformation (statistics and dynamics) 
PACS 64 . 70 . Nd - Structural transitions in nanoscale materials 
PACS 87 . 15 . ad - Analytical theories 

Abstract. - We study the structures of a ribbon or ladder polymer immersed in poor solvents. The 
anisotropic bending rigidity coupled with the surface tension leads ribbon polymers to spontaneous 
formation of highly anisotropic condensates in poor solvents. Unlike ordinary flexible polymers 
these condensates undergo a number of distinct layering transitions as a function of chain length 
or solvent quality, and the size of condensates becomes non-monotonic function of chain length. 
We show that the fluctuations of the condensates are in general small and these condensates are 
stable. 
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Recently considerable attention has been devoted to the 
physics of semiflcxible polymers, i.e., chains which have 
significant bending stiffness. These kinds of polymers ex- 
hibit interesting physics in the form of liquid crystalline 
phases and are of great importance in biology and bio- 
physics. For example, DNA is the most well-known semi- 
flexible polymer. Much of the interest has focused on 
semiflcxible chains with isotropic bending elasticity [1—7]. 
Macroscopically this corresponds to modeling the chain 
as a cylinder with a circular cross section. In many cases, 
however, the cross section is highly anisotropic [8-12]. 

The anisotropic polymer can be modeled as a ribbon, of 
length L, width w, and thickness t assuming L 3> w > t 
(fig. 1) with anisotropic bending elasticity. Indeed, if the 
ribbon is composed of an isotropic material, the ratio of 
the elastic constants for bend in the easy and hard direc- 
tions is e e /eh — (t/w) 2 [13] . This strong dependence on 
the ratio of t/w implies that the ribbon polymer shows a 
very anisotropic elastic response. On the chemical scale, 
the chains do not normally consist of an isotropic material, 
and details of chemical bonding are important for bend, 
thus, even higher ratios of e e /et, are possible. 

Some of the solution properties, in particular the liquid- 
crystalline behavior, have been studied theoretically [10]. 
Here we study the structure of a single ribbon polymer im- 
mersed in a poor solvent. Assuming the system is well be- 




Fig. 1: Model for anisotropic semiflexible polymer 



low the O temperature, the chain collapses completely and 
forms a condensate which is effectively a polymer melt. 

For isotropic semiflexible chains such as DNA, the con- 
densate often forms a toroid [1,2, 14-21] and sometimes a 
globule [22-25] through intermediate states [28-30]. For 
ribbon polymers, we show in this paper that the conden- 
sates forms an annulus. The annulus formation manages 
to avoid bend in the hard direction, while reducing con- 
tact of the chain with the solvent. The novel character 
of these annular condensates is that multilayer annuli can 
form and they undergo sudden changes as a function of 
chain length or surface tension. This is in marked con- 
trast to the case of ordinary flexible or semiflexible poly- 
mers where no dramatic changes occur. This quantization 
is a direct consequence of the anisotropic elasticity. 

To model the system, we introduce two surface tensions 
and 7 e which corresponds to the surface energies of con- 
tact between the top (narrow surface) of the polymer and 
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Fig. 2: Annular condensate for short ribbon polymer in poor 
solvent 

the solvent, and the surface energies of contact between 
the side (fiat and wide surface) of the polymer and the 
solvent, respectively. If the chain consists of anisotropic 
molecules, any ratio of ^e/lh is possible in principle. In 
this paper, however, we assume 7 e /7/i ~ 1- 

A large number of parameters appear in this model. 
Once final results have been obtained, it is useful to ignore 
numerical prcfactors and to make some crude estimates. 
Then, we write w ~ t w a and 7 e ~ ~ gkT/a 2 , where g 
is of order unity and a is a length of order of one Angstrom. 
It is convenient to define two bare persistence lengths, l e = 
e e /kT and lh = Sh/kT which are essentially the scales 
on which a free chain bends in the two directions due to 
thermal fluctuations. 

We begin by examining the shape of condensates when 
we slowly lengthen the chain. For short chains, bend costs 
clastic energy, then the chain more-or-less remains as a 
rod. For slightly longer chains, however, the overlap in- 
duced by bend compensates the elastic energy by reducing 
its surface energy, therefore, the chain may form a ring. 
For a circular ring of radius R it is easy to show that the 
bending energy is -^e e L/R 2 . The change in surface energy 
in bending into a ring is AnRw^/e — 2Lwj e . Minimising 
it over R and equating the free energies gives the critical 
length, L* = (3%/3tt/2) y/e e /{-y e w) at which ring forma- 
tion occurs. The chain overlaps by 1/2 turn, at L = L*. 
Our interest here is in the case L ^ L* so the chain winds 
around many times. 

The simplest way the chain can pack is to form a disk 
with a hole in the middle (annulus) as shown in fig. 2. This 
annulus has height w and inner and outer radii Ri and 
i? , respectively. The total energy consists of three terms: 
(1) surface energy of the top and bottom surfaces, (2) 
surface energy of the inner and outer exposed sections of 
the annulus, (3) bending energy of the chain. This last 
term is due to bending only in the easy direction. 

Because the volume of the chain is conserved, R and 
Ri are related by 

tLw = wk(R 2 - R 2 ). (1) 

There is thus only one variable to minimize over (say Ri). 
Also, since the area of the upper and lower surfaces is 
fixed (i.e., both are always of area tL), there are really 
only two relevant terms in the total energy. The lateral 



surface energy is easily calculated as 

2ir le w(R + R l ). (2) 
The bending energy is calculated by 

where s is the arc length and 6 is the angle between the 
tangent of the ribbon curve and some fixed direction in 
the plane of the disk. 

In one complete turn of the ribbon, 9 increases by 27r 
and s increases by 2-7rr where r is the distance of the chain 
from the center of the annulus. This gives dd/ds = 1/r. 
During the same change in s the radius increases by one 
ribbon thickness, i.e., dr/ds = t/(2Ttr). This allows us to 
write the bending energy as 7re e t -1 \n(R /Ri). The total 
energy is then 

U = "^(if) + 2 ^M R o + R l ), (4) 

subject to R 2 = R 2 + tL/%. 

It is convenient at this point to introduce two charac- 
teristic lengths. One is A = e e /(2j e wt), and the other is 
r = y/ Lt/ir. For long chains, T is much larger than A, and 
we shall see Ri « A and R ~ T. In this limit, A <C T, 
the annulus is "fat", i.e., it has a small hole in the middle. 
In the opposite limit, A ^ T, the hole is large and the 
annulus is "thin" . It looks like a ring. 

In the long-chain limit, Ri <C T, the total energy is 
approximated by 

7T£ e / ItRi \ 

U = const. + 2ir~/ e wRi In I —= I , (5) 

t \VLtJ 

which has a minimum at R* = A. Using our crude 
estimates give R* « l e /g, i.e., the internal radius is 
somewhat less than the bare persistence length. In this 
limit the internal radius of the annulus is independent of 
chain length, whereas the external radius is approximately 
R ~ (tL/n) 1 / 2 . A more accurate expression for R* can 
be found by expanding in ascending powers of a = A/T. 
This gives R* = A(l — a + |a 3 — ...). From this it is 
clear that, even though R* becomes asymptotically A, the 
approach is rather slow as L~ x l 2 . 

In the limit A » T, we can again expand the total 
energy and find the minimum at R* = (r 2 A/2) 1 / 3 , or 
more crudely, R* w {Ll e a/g)^ 3 . The crossover between 
the thin and fat regimes occurs at T = A. A crude estimate 
of this gives L c = a(L*/a) 4 . This implies that the chain 
must be rather long to be in the fat regime. 

The results just derived correspond to the case of one 
layer. The novel point about this system is that it can 
choose to form more than one layer. This decreases the 
surface area exposed by the top and bottom surfaces but 
comes at the expense of extra surface energy on the inside 
and outside surfaces and extra bend both in the hard and 




Fig. 3: tubular condensate for long ribbon polymer in poor 
solvent 



easy directions. For second layer to form the chain must 
bend in the hard direction. The best way of doing this is 
to join the two layers together at the outer radius (fig. 3). 

We need to calculate the bending energy of the chain in 
the hard direction. The chain must bend between the two 
layers. It must jump a distance w in the vertical direction 
while winding a distance Q = 2irR in arc length. If we 
let cj) be the angle made by the chain with the plane of the 
lower annulus then the bending energy can be written as 



(6) 



This is to be minimized subject to a constraint that the 
chain jumps up a distance w, i.e., J Q ds(f> = w. Here we 
have made the very good approximation that the angle <p 
is always small. Using the method of Lagrange multipliers 
and the Eulcr-Lagrangc equations we find that the mini- 
mum energy trajectory coresponds to the solution of the 
differential equation: d 2 (f>/ds 2 — const. Using the bound- 
ary conditions (f)(0) = = 0(0) the constraint gives 



and the bending energy becomes 

6e h w 2 



U h = 



O 3 



(7) 



(8) 



It is fairly clear from this formula that the chain wants to 
spread the bend over as much length as possible, so that 
it is best to bend at the outer radius of the annulus and 
Q = 2itR . 

We now examine the transition from one to two layers. 
The one-layer energy is 

E/i = ?Ef£ hi f^L \ + 2ir Je w(R + R t ) + 2 7h Lt, (9) 



subject to R 2 = R 2 + tL/n. The two- layer energy is 



U 2 
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(10) 



subject to R 2 = R 2 + tL/2n. Here jhLt is the surface 
tension of the "top" of annulus. 

The transition from one to two layers occurs when U\ = 
Ui- To solve this equation we need to minimize U\ over 
Ri and similarly for XJi- In principle different Ri will be 
associated with the two geometries. However, in the limit 
where Ri <C (Lt) 1 / 2 the result is the same for both, i.e., 
R 2 ps e e /{2-f e wt). This gives 
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where n is the number of layers (1 or 2). Equating the two 
energies for n = 1 and n = 2, and ignoring the logarithmic 
factors, gives 



-j h Lt + 2{V2- l)7rw 7e 



3e h w 2 
~V2n 3 



-3/2 



(12) 



The four terms in this equation represent the four kinds 
of energy differences between the two layer and one layer 
system. In order these are: Easy bend; top/bottom sur- 
faces; side surfaces; hard bend. The driving force for the 
transition to two layers is the top/bottom surface energy. 
This is opposed by the other three terms. We are mainly 
interested here in the case where the hard bend is the 
dominant penalty. In this case the transition occurs when 



L' 



(13) 



where we have ignored numerical prcfactors. Even more 
approximately this can be written L ps a(lh/a) 2 ^ 5 where 
lh is the persistence length in the hard direction. Provided 
lh S> l e (it is true when L 3> R4) our approximations are 
valid. 

At the transition the jump in Ri is small, since Ri is 
about the same for both layers. However, there is a sudden 
jump in R by a factor of l/v2- In the case of a three 
layer system, the chain must jump once on the outside 
radius and once on the inside radius. Since the jump on 
the inside radius must take place over a shorter amount of 
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Fig. 4: Layering transition in a fat condition F/A ~ 10. Fhe 
horizontal axis is the polymer length L/w. The upper and the 
lower curves represent R and Ri respectively and the middle 
step line represents the height of the tubular condensate in unit 
of ribbon width w. 
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Fig. 5: Layering transition in a thin condition F/A ~ 0.5. The 
horizontal axis is the polymer length L/w. The upper two 
curves represent R and Ri respectively and the lower step 
line represents the height of the tubular condensate in unit of 
ribbon width w. 



arc length, its energy cost is expensive. The energy for n 
layers is then 



27r 7 e J 2nRl ds\(j)(s)\r(s) . The total energy penalty is then 



U[r 



nne e ( R Q 



+n(27T7 e w)(i? +Ri 



(2tt) 3 



R„ 



2 T 

-JhLt 

n 

n-1 



(14) 



with R = (R 2 +n 1 tL /it) 1 / 2 , and where [x] indicates the 
integer part of x. 

In order to show the layering transition of equilibrium 
form in a fat condition, in fig. 4 we plot the outer radius 
of the condesates R , the radius of the inner hole of the 
condensates Ri, and the quantized height of the conden- 
sate nw. This fat condition is always achieved in the long 
chain limit. For large L, the size of the condensate is 
roughly proportional to \f~L. Below the medium chain 
length, there exists a certain parameter region where the 
condensate is thin. We show in fig. 5 behavior of the lay- 
ering transition in a thin condition. 

We have assumed here that the system can form a well- 
defined annulus. It is important and interesting to exam- 
ine the fluctuations around the equilibrium shape of the 
system. We will concentrate here on the case of a single 
layer. There are two types of fluctuations. The first kind 
is fluctuations of the shape in the radial direction, i.e., 
distortion from a circle. We discuss it later. 

The second kind is fluctuations out of the plane which 
cause two kinds of energy penalty. The first is arises from 
hard bend and the energy penalty is e^/2 ds(d<f)/ ds) 2 , 
where (f> is the angle made by the tangent to the rib- 
bon out of the plane of the annulus. The second arises 
from the surface area exposed due to differences in the 
height z(s) of the layers. This gives rise to an energy 
penalty 7 e ds\z(s) — z{s — 27rr(s)}|. Noting that 
z(s) = Jg ds'<fi(s') and expanding for small gradients gives 



■ 27T7 e 



ds r(s)\<j>{s)\. (15) 



This kind of total energy is peculiar because of the form 
of the surface tension term, |</>(s)|. This term is first order 
and not differentiable about the equilibrium position <fi = 
0. This means that the usual eigenmode techniques we 
might have used to analyse the system are not applicable. 

However, a good estimate of the size of the fluctuations 
and in particular 5z = y/ ((z(L) — z(0)) 2 ) can be found 
by mapping the problem onto that of a semiflcxiblc chain 
in a nematic field. From now on we ignore all numeri- 
cal prefactors. For our purposes it is convenient to use 
an approach developed by Odijk [31,32]. He considers 
a wormlike chain directed by a parabolic potential /x</> 2 , 
where [i is a constant. The chain has a typical angular 
deviation 5(f) = y (tfi 2 ) . Such a chain can be viewed as 
a series of steps of size Am = lh(5(f>) 2 in the plane of the 
annulus and X z = Am (50 = lh(5<fi) 3 perpendicular to the 
annulus. A chain of length L undergoes an anisotropic 
random walk of L/Xu steps in the plane, with each step 
travelling a distance X z perpendicular to the plane. This 
implies the deviation out of the plane is Sz = \J Llh(S(f>) 2 . 

Odijk shows that the strength of the potential and the 
angular deviations are related by (S(f>) 4 = kT/{fxlh), so 
that Sz = y/kTL/ 1 [i. This result cannot be applied di- 
rectly to our problem, since we have a confining potential 
of the form r y e r(s)\(f>(s)\. However, the details of the shape 
of the potential cannot matter for the gross behaviour of 
the system. In particular we can replace our potential by a 
harmonic one, provided we choose p, self-consistcntly, i.e., 
we choose so that for a typical value of the angular de- 
viation the two potentials coincide. Equating the two po- 
tentials gives 7 e r(s)|0| = /i</> 2 . This together, with the re- 
lation (<50) 4 = KT/(jjl h ) gives y, = ( 7e r( S )) 4 / 3 (/ ft /fcT) 1 /3. 
Note that since r varies along the chain, the effective con- 
fining potential /i also varies along the chain. We can now 
integrate over the whole annulus to get the variation in 
height Sz. Between r and r + dr there is a total length 
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rdr/t and hence the square of the step in the z direction 
is dz 2 = (kT '/ 1 {'-y e r)) i / d 'l h 1 t' i t~ 1 rdr. Integrating over the 
whole annulus gives: 

JC kT h /„ 2/3 D 2/3 

Note that this deviation grows only very slowly with the 
chain length. In, particular, for very long chains we have 
Sz = (fcT/( 7e ? /l )) 2 / 3 (/ /l /0 1/2 t 1/6 i 1/6 . Setting t = a and 
7 e = gkT/a 2 gives Sz « ag^ 2 ^{L jlh) 1 ^ . Thus, in general 
the overall fluctuations will be less than, or of the same 
order as the annular thickness, w, and the annulus will be 
well-defined and relatively flat. This is the conclusion for 
one layer. 

For many layers this conclusion remains true, since in 
that case the effective persistence length scales with the 
number of layers and the radius also scales with the num- 
ber of layers. This implies that the absolute deviation in 
thickness remains roughly constant and the relative vari- 
ation decreases as more layers are added. 

We now discuss fluctuations of radial direction. The 
simplest mode is one where the inner radius changes from 
its equilibrium value, R*. In the long-chain limit, where 
R* <C vLt the energy penalty for deviating from equilib- 
rium is 

w _~(*zsy. (17) 

Equating this to kT/2 gives a typical size of fluctuation 
y/({Ri - R*) 2 ) = Ay/t/(Trl e ). For l e /t = 10 this gives a 
relative fluctuation of 13%. We expect that the fluctua- 
tions within the plane will be reasonably small. 

In conclusion, ribbon polymers which have large bend- 
ing rigidities in lateral direction spontaneously form 
anisotropic condensates in poor solvents. The size of those 
peculiar condensates quantized as a variation of the length 
of polymers and solvent quality. 

When the material is homogeneous, the surface tensions 
and elastic energies are determined by geometry of the 
ribbon. In general, for homogeneous materials, we may 
assume e e < eh and ~f e « jh, which is the case we dis- 
cussed in this paper. On the chemical scale, however, 
other anisotropic conditions are possible in principle. The 
ribbon polymers with e e < and wj e < tjh [33], form a 
different type of anisotropic condensates in poor solvents, 
that will be discussed somewhere else. 

* * * 

The authors thank Jean-Louis Sikorav for discussions 
and comments on the manuscript. This work was initi- 
ated with support by NTT Basic Research Laboratories. 
Y.Y.S. wishes to acknowledge CEA, IPhT and Takushoku 
University, RISE for financial support. 



REFERENCES 

[1] Klimenko S. M., Tikkchonenko T. I. and Andreev 

V. M., J. Mol. Biol, 23 (1967) 523 
[2] Grosberg A. Yu., Biophysics, 24 (1979) 30 
[3] Frank-Kamenetskii M. D., Anshelevich V. V. and 

Lukashin A. V., Sou. Phys. Usp., 30 (1987) 317 
[4] Marko J. F. and SlGGlA E. D., Macromolecules, 27 

(1994) 981 

[5] Odijk T., J. Chem. Phys., 105 (1996) 1270 
[6] Marko J. F., Phys. Rev. E, 55 (1997) 1758 
[7] Haijun Z., Yang Z. and Zhong-CAN O.-Y., Phys. Rev. 

Lett, 82 (1999) 4560 
[8] Schluter A.-D., Adv. Mater., 3 (1991) 282 
[9] Nyrkova I. A., Semenov A. N., Joanny J.-F. and 
Khokhlov A. R., J. Phys. II (Paris), 6 (1996) 1411 
[10] Nyrkova I. A., Semenov A. N. and Joanny J.-F., J. 

Phys. II (Paris), 7 (1997) 825 
[11] Nyrkova I. A., Semenov A. N. and Joanny J.-F., J. 

Phys. II (Paris), 7 (1997) 847 
[12] Golestanian R. and Liverpool T. B., Phys. Rev. E, 
62 (2000) 5488 

[13] Landau L. D. and Lifshitz E. M., Theory of Elasticity 

(Pergamon, London) 1970, p. 77 
[14] Bloomfield V. A., Biopolymers, 31 (1991) 1471 
[15] Bloomfield V. A., Biopolymers, 44 (1997) 269 
[16] Li A.-Z., Fan T-Y. and Ding M., Science in China B, 

35 (1992) 169 

[17] Ubbink J. and Odijk T., Biophys. J., 68 (1995) 54 
[18] Hud T., Downing K. H. and Balhorn R., Proc. Natl. 

Acad. Sci. USA, 92 (1995) 3581 
[19] Park S. Y., Harries D. and Gelbart W. M., Biophys. 

J., 75 (1998) 714 
[20] Pereira G. G. and Williams D. R. M., Biophys. J., 80 

(2001) 161 

[21] Takenaka Y., Yoshikawa K., Yoshikawa Y., 
Koyama Y. and Kanbe T., J. Chem. Phys., 123 (2005) 
014902 

[22] Grosberg A. Yu. and Zhestkov A. V., J. Biomol. 

Struct. Dyn., 3 (1986) 859 
[23] Vasilevskaya V. V., Khokhlov A. R., Kidoaki S. and 

Yoshikawa K., Biopolymers, 41 (1997) 51 
[24] Jary, D. and Sikorav J.-L., Biochem., 38 (1999) 3223 
[25] Odijk T., Phil. Trans. R. Soc. Bond. A, 362 (2004) 1497 
[26] Cooke I. R. and Williams D. R. M., Physica A, 339 

(2004) 45 

[27] Yoshikawa Y., Yoshikawa K. and Kanbe T., Biophys. 

Chem., 61 (1996) 93 
[28] Schnurr B., Mackintosh F. C. and Williams D. R. 

M., Europhys. Lett, 51 (2000) 279 
[29] Pereira G. G. and Williams D. R. M., ANZIAM J., 

45 (2004) C163 
[30] Yoshinaga N., Yoshikawa K. and Ohta T., Euro. 

Phys. J. E, 17 (2005) 485 
[31] Odijk T., Macromolecules, 16 (1983) 1340 
[32] Odijk T., Macromolecules, 19 (1986) 2313 
[33] Nelson J. O, Saven J. G., Moore J. S. and Wolynes 

P. G., Science, 277 (1997) 1793 



